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tential of the Debye- Yukawa form are considered. The pressure is ob- 
tained from the stress tensor of the field corresponding to the Yukawa 
interaction, by a suitable statistical average. This approach is espe- 
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1. INTRODUCTION 

In a recent paper [|l| it has been discussed how the pressure in a Coulomb fluid can be 
obtained by a suitable statistical average of the Maxwell stress tensor. The present 
paper is an extension of these calculations to the case of a fluid made of particles with 
Yukawa interactions. These Yukawa fluids have been recently reviewed 0,^]. They 
are of interest because they are often used as simplifled models for some complex 
Coulomb fluids (for instance dusty plasmas): one species S of charged particles is 
given a microscopic description, while the other charged particles are taken into 
account only through the screening that they cause to the interaction between the 
^'-particles. One is left with one species interacting through a screened Coulomb 
potential which is assumed to be of the Debye- Yukawa type. In the following, this 
model will be considered for its own sake. In particular, the case when the fluid is 
in a curved space (of interest both for numerical simulations and for looking at the 
curvature effects) will be investigated. 

To start with, the system under consideration is, inlR^, a gas of point-particles of 
number-density n, each of charge q, with a pair interaction v{r) = q^G{r) depending 
on the distance r, where G{r) has the Yukawa form 

G(r) = ^"P(-"^^ (1.1) 

r 

which is the Green function of the Helmholz equation. From the thermodynamic 
expression for the pressure 

P = -dF/dV (1.2) 

where F is the free energy and V the volume, one can derive |^ the standard virial 
expression for the excess (over ideal) pressure: 

"ex = --n q / r—g[r)dv (1.3) 

where g{r) is the pair distribution function (in the present paper, integrals without 
domain speciflcation are meant to be extended to the whole manifold). In Section 2, 
it will be shown how (1.3) can be alternatively derived from the stress tensor of the 
field which carries the Yukawa interaction. In Section 3, the above considerations 
will be extended to the case of a system confined on the surface S^ of an hypersphere. 



Section 4 deals with the two-dimensional analogs of these systems. The limiting case 
of Coulomb systems is considered in Section 5; the case of a pseudosphere (a surface 
of constant negative curvature) is revisited in Section 6. The results are discussed 
in Section 7. 
2. YUKAWA GAS INIR* 

In terms of the Yukawa field (f){r) created by some charge distribution, the energy 
density is (l/87r)[(V0)^ + a;^0^] and the corresponding stress tensor is P] 

T,. = ^[d,<Pd,<P - ^5^.((V0)2 + aV')] (2.1) 

where the Greek indices label the three Cartesian axes x, y, z. In terms of the particle 
positions r,, 

0(r) = g^G(|r-r,|) (2.2) 

i 

The excess pressure is the negative of the statistical average of any diagonal 
element, say — < Tj.^. >, at some point which can be chosen as the origin. Following 
the same steps as in ref. [|I|, one can decompose the excess pressure into nonself and 
self parts. With the rotational symmetry taken into account, the nonself part can 
be written as 

Pnonseif = '^ j dv^dv^[VG{n) ■ VGir^) + 3a^G{n)G{r2M\r2 - r,\) (2.3) 

Writing each G as the Fourier transform of 47r/(A;^ + a^), using as integration vari- 
ables ri and r = r2 — ri, and performing first the integration on ri , one recovers 
(1.3). As to the self part of — < T^^. >, 

P.eif = -^ / dr[{dMr)f - {dyG{r)f - {d,G{r)f - a' {{G{r)f] (2.4) 

it is a divergent integral (at small r), which can however be regularized through the 
use of the same physical argument as in ref. 0] : for computing in an unambiguous 
way the force per unit area across the yOz plane, one must assume that no particle 
sits on that plane. Thus, one removes from the integration domain in (2.4) a thin 
slab —£< X < e and takes the limit e — * at the end of the calculation. After some 
algebra with Fourier transforms, this prescription gives Pseif = 0. 



Thus , inE^, 

-'ex -* nonself V-^''-*/ 

and the stress tensor approach simply reproduces the standard virial formula (1.3). 
3. YUKAWA GAS IN S=^ 

We now consider a Yukawa fluid living on the three-dimensional "surface" 5''^ of 
an hypersphere of radius R. The Yukawa potential (1.1) now must be changed into 
the Green function of the Helmholtz equation on 5''^ which is 

^, , 1 sinhuj (it — ih) ^ 

G(iP) = ^ , for aR > 1 3.1a 

R sm ip smh un 

^, ,, 1 sint^(7r - V^) 

<^(^) = t; ^ fo^^ aR< 1 (3.16) 

R sin ^ sin tun 

where ip is the angular distance seen from the center of the hypersphere [Rtp is the 
geodesic distance) and cu = \a'^R'^ — 1|^/^. 
3.1. Pressure from the stress tensor 
The analog of (2.3) is 

Pnonself = ^J dVidV2[V qG (^oi) " VoG(^02) + 3a2G'(^oi)G'(^02)]^(^12) (3.2) 

where is the arbitrary point at which Pnonseif is evaluated, dVi is a volume element 
around point i, and Ripij is the geodesic distance between points i and j. After a 
calculation described in Appendix A, (3.2) becomes 

Pnonseif = " ^^'^ '/ ^^^^^"^^"^^ ^^'^^ 

where G{ip) must actually be regarded as a function of two independent variables: 
the angular distance ip and the hypersphere radius R, as seen on (3.1). In (3.3), the 
volume element is dV = AttR^ siii^ ipdip. This result (3.3) is the generalization of 
(1.3) to a curved space (the hypersphere). It should be noted that now the distance 
r has been replaced by the radius of curvature R. Only in the flat system limit does 
G{ip, R) become a function of the sole variable r = Rip and RdG/dR = rdG/dr. 

Furthermore, now, the stress tensor approach provides another contribution to 
Pex'- the properly regularized part Pgeif does not vanish in the S^ case. Indeed, the 
analog of the integral (2.4) can be split into two pieces Pq and Pi corresponding to 



the contributions of geodesic distances to the origin smaller and larger than Ripo, 
respectively. Only Pq is divergent (at the origin) and must be regularized by the 
same prescription as in the IR^ case. It is convenient to choose an infinitely small 
i/jq. Then, the regularized Pq can be computed by using the small-?/' form of (3.1), 
which is just the Coulomb potential inlR^. This calculation has been done in ref. [|T| 
with the result 
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nq 







6Ripo 



+ 0(V^, 
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(3.4) 



As to Pi, with the rotational symmetry taken into account, it can be written as 

2 



Pi 



nq 



1 dG' 



3a'G 



2^2 



dV 
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This integral (3.5) is computed in Appendix A, and the final result is 



(3.5) 



self 



Po + Pi 



nq 



TTO^R cot TTUO 



sm TTUJ 



Rtu 



for aR < 1 



(3.6) 



For aR > 1, u must be replaced by iu in (3.6). 

Thus the stress tensor approach gives a total excess pressure Pex which is the 
sum of the nonself part (3.3) and the self part (3.6). While the nonself part is the 
analog of the total excess pressure (1.3) in IF^, there is now in S^ an additional 
non- vanishing self term (3.6). 
3. 2. Pressure from the free energy 

Another way of defining and computing the pressure in the case of 5*^ is to use 
the standard equation (1.2). If the total potential energy is 



i<j 



(3.7) 



the excess free energy is given by 



/3Pex = -Inl— I dVi... dVN exp 



-(3Y.G{^^ 



u; 



i<j 



(3.8) 



where /3 is the inverse temperature and A^ the number of particles. From (1.2), 
where here V = 27r^P^, one finds for the excess pressure the value (3.3) which was 
the nonself part in the stress tensor approach. 



The self part (3.6) appears only if one includes in the total potential energy 
the self-energies of the particles. Up to some (infinite) volume-independent additive 
constant, the self-energy of a particle is obtained in Appendix A as 

Cseif = ~7^(f'~B cotvrcj for aR < 1 (3.9) 

2 R 

If one adds to the free energy the self term A^Cgcif, one does obtain from (1.2) the 
additional self contribution (3.6) to the pressure. 

These results will be discussed in Section 7. 
4.TWO-DIMENSIONAL YUKAWA FLUIDS 

In this Section, the above considerations are adapted to the case of two- 
dimensional systems. These two-dimensional systems are toy models of theoretical 
interest. It should be remembered that, in many experimental situations of charged 
particles confined on a surface, the interaction nevertheless is the three-dimensional 
Coulomb potential rather than the two-dimensional interactions which are consid- 
ered here. 
4. I.Yukawa gas inE^ 

InE^, the Green function of the Helmholz equation is 

G{r)=K^{ar) (4.1) 

where Kq is a modified Bessel function and (1.3) is replaced by 

Pex = --n\^ J r—g{r)dr (4.2) 

In the field approach, the energy density is (l/47r)[(V0)^ + q;^0^] and the corre- 
sponding stress tensor is 

T,. = ^[d,<Pd,<P - ^V((V0)' + a'<P')] (4.3) 

Zn Z 

where the Greek indices now label two Cartesian axes x, y. 

When the rotational symmetry is taken into account, the analog of (2.3) has the 
simpler form 

2 2 

Pnonscif = ^a^ I dr^dr2G{n)Gir2)gi\r2 - r^\) (4.4) 



(in two dimensions, the derivatives of G cancel out). As in the three-dimensional 
case, the expression (4.4) of Pnonseif can be brought to the form (4.2) while the prop- 
erly regularized Pseif vanishes. Thus, Pcx = -Pnonseif, and the stress tensor approach 
reproduces (4.2). 
4. 2. Yukawa gas in S^ 

As shown in Appendix B, the Green function of the Helmholz equation on a 
sphere S"^ of radius R is 

G{e) = - ^ ^ P,{-cose) (4.5) 

2 sm UTT 

where 9 is the angular distance, P^ is a Legendre function, and 
u = (1/2) [-1 + (1 - 4a2p2)i/2] (^if 2aR > 1, z/ is complex). 
The analog of (4.4) is 

Pnonself = ^^^ J dS,dS2G{9o,)G{9o2)g{ei2) (4.6) 

where dSi is an area element around point i. After a calculation similar to the one 
for S3 in Appendix A, (4.6) becomes the analog of (3.3) 



nonself 



-\nV J R^gi9)dS (4.7) 



{G{9) depends also on R through u). When Pseif is split into the contributions Pq 
(Pi) of geodesic distances smaller (larger) than ^6*0, both parts remain finite as 
^0 -^ 0. The regularized Pq is the same as for a plane two-dimensional Coulomb 
system [|l| (Pq = — ng^/4) and 

2 

P^ = !^«2 hc^o^fds (4.8) 

Atx J 

The integral in (4.8) involves ^_^[Pp{x)Y'dx which is tabulated in ref. P], and one 
obtains 
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Po^Pi 



4 2u + l 



7r2 



sin^ TTZ/ 



2^'(z/ + l) 



(4.9) 



where ip' is the derivative of the psi function (the psi function being the logarithmic 
derivative of the gamma function). 



The pressure can also be derived from the thermodynamic relation (1.2) (with the 
volume V replaced by the sphere area S). Again, if only the two-body interactions 
(f'GiOij) are used for defining the free energy F, only Pnonseif is obtained, and, for 
Pseif to appear, it is necessary to add to the free energy the self-energies of the 
particles. Each particle is found (see Appendix B) to have, up to some (infinite) 
volume-independent additive constant, the self-energy 

(4.10) 



_ ^ 

Cself — — 



TT 
In i? — ^(z/ + 1) cot TTZ/ 



This gives to the pressure the self contribution (4.9). 
4. 3. Yukawa gas in a pseudosphere 

The pseudosphere is a surface of constant negative curvature. Since it is infinite, 
on it it is possible to study systems which are both infinite and curved. 

Let a be the "radius" of the pseudosphere, such that the Gaussian curvature is 
— 1/a? (instead of 1/ R^ on a sphere) and let ar be the geodesic distance (instead 
of R9 on a sphere). As shown in Appendix C, the Green function of the Helmholtz 
equation now is 

G(r) = g,(coshr) (4.11) 

where Qu is a Legendre function of the second kind and 

v= (l/2)[-l + (l + 4aV)i/2]. 

The nonself pressure is given by (4.6) where G now is G{t) as given in (4.11). 
After a calculation similar to the one in Appendix A, one finds the analog of (4.7): 

Pnonseif = -\n\^ j a^^^g{T)dS (4.12) 

As to the self pressure, its part Pq is unchanged (Pq = —nq^ /A) while Pi is given 
by (4.8) where G is G{t). The integral involves J^[Qu{x)]'^dx which is tabulated in 
ref. [HI, and one obtains 



Pscif — Pq + Pi . 



"^'^l + 2aV^'^" + '^ 



2i/-M 



(4.13) 



5. COULOMB LIMIT 

The results of ref. [jl| for a Coulomb fluid, the one-component plasma, can be 
retrieved in the limit a ^ of the Yukawa fluid. However, before taking that limit. 



one must add a neutralizing background to the Yukawa fluid. The presence of a back- 
ground is taken into account by changing g into h = g — 1 in (1.3), (3. 3), (4. 2), (4. 7). 
5.1. Coulomb limit inE^ 

As a ^ 0, G ^ 1/r and (1.3) becomes 

1 /■ 1 

^ex = ^nonseif = -n^Q^ / -h{r)dr (5.1) 

J r 

i.e. the excess pressure is one third of the potential energy density. 
5. 2. Coulomb limit in S^ 

As a — * 0, 

Since the system is finite, 

n f h{ip)dV = -1 (5.3) 

and (3.3) becomes 

Pnonscif = l:nV I \[{ir- ^) cot^ - \]h{i,)dV - ^^ + o(l) (5.4) 

Furthermore, as a — ^ 0, (3.6) becomes 

Thus, one retrieves the result of ref. [^ for the total excess pressure in the Coulomb 
case a = 0: 

1 /" 1 1 no 

Peyi = ^nonself + ^self = -n^Q^ / —f^li'^- ^) COt V' " -]h{ij)dV - — — (5.6) 

It should be noted that, in the limit a — > 0, Pnonseif and Pseif have opposite 
divergent terms 0(a~^) which cancel each other in their sum Pex- Thus, it is essential 
to keep the self term for retrieving the finite result (5.6). 
5. 3. Coulomb limit inH^ 

As a ^ 0, KQ^ar) = — ln(ar/2) — 7 + o(l) (where 7 is Euler's constant) ||^ and 
rdG/dr ^ — 1. Since the Coulomb fluid exhibits perfect internal screening, i.e. 

n / h{r)dv = —1 (5.7) 



(4.2) gives the known explicit exact result [0 



3 - P -_!!^ r^^ 

ex -'■ nonself . yo.Q 



5. 3. Coulomb limit in S^ 

As a ^0, z/ = -a^R^ + o{a^), P^{-cos9) = 1 + 2z/lnsin(^/2) + o(a2) [|, and 
RdG/dR = -{l/a^R^) + o(l). Since the system is finite, 



n / h{e)dS = -1 (5.9) 



and (4.7) becomes 
while (4.9) becomes 



4a2i?2 
nq 
Thus, one retrieves the result of ref . M : 



nq'^ 

nonsolf = - , or^o + o{l) (5.10) 



Pself - ^ 



2 r 1 



^ + ^+^(^) 



(5.11) 



+ P.cif = -^ (5.12) 

Here too, Pnonseif and Pgeif have opposite divergent terms Oi^a'"^) which cancel 
each other. 
5. 4. Coulomb limit in a pseudosphere 

As a ^ 0, z^ -^ 0, G{t) —>■ (5o(coshr) = — lntanh(r/2), and adG/da -^ 0. Thus 
-Pnonseif = 0. On the other hand, (4.13) becomes Pgeif = —nq^ /A. Thus, 

2 
-Pex = -Pself = -7- (5.13) 

Here, the excess pressure entirely comes from the self part. 
6. VIRIAL EXPANSION IN A PSEUDOSPHERE 

In a pseudosphere, Jancovici and Tellez have defined an excess pressure by 

the virial expansion in powers of the density n 

00 

/3pex = E^fc^' (6-1) 

fc=2 

Obviously, (6.1) does not agree with (5.13). The present Section is an erratum to 
ref. 0: (6.1) is not a good definition of the pressure. The considerations about the 
virial expansion in ref. should be replaced by what follows. 

10 



Since on a pseudosphere, as the size of a domain is increased, its perimeter 
grows as fast as its area S, there is no unique thermodynamic hmit for the free 
energy per unit area F/S. A reasonable definition of a bulk quantity is provided 
by the usual virial expansion of the free energy, obtained by manipulations of the 
partition function: the nonself part of the excess free energy per particle can be 
defined by 

CO 75 

/3/nonself = E IT^^'"' (6-2) 

k=2 '^ ^ 

where each virial coefficient B^ is to be computed in the infinite system limit (for 
eliminating the boundary effects) before the sum is performed. 

From this free energy (6.2), one can derive the nonself part of the excess pressure 
by the standard relation Pnonseif = n'^dfnonseXildn provided one takes into account an 
unusual feature (which has been disregarded in ref. 0): the interaction law G and 
therefore the excess free energy per particle /ex depends on the radius of curvature 
a (the virial coefficients B^ do depend on a): it is convenient to consider /ex as a 
function of the two independent variables n and no? . The variation of density in the 
above definition of the excess pressure can be considered as associated to a variation 
of a while the average number of particles in some domain remains constant as the 
area of this domain varies proportional to a?. In other words, for defining the excess 
pressure, the partial derivative of the excess free energy with respect to n must be 
taken at constant na^, and one must write more precisely 

"nonself = U I I (6.3) 

\ / na^ 

In the case of a Coulomb system of point particles, for dimensional reasons, /nonscif 
depends on n and a only through the combination no? {Bk is proportional to a^^'^"^-'), 
and (6.3) gives Pnonseif = 0. 

Furthermore, self effects must be taken into account. For a Coulomb system, 
with an interaction G = — lntanh(r/2), the self-energy of a particle of small radius 
ro = aro 

eseH = /se,f = -yln- = --ln^-, (6.4) 



11 



does give a contribution to the pressure 

dn I n 4 



U^n'i^^] =-'4 (6-5) 



in agreement with (5.13). 
7. DISCUSSION 

In a flat space ]R^ or I?, the stress tensor approach simply reproduces the expected 
pressure. There is no self contribution. 

The situation is more involved in the finite curved spaces S^ and S'^, and in 
the infinite pseudosphere. The very definition of the pressure is not obvious. An 
operational definition would be the force per unit area exerted on a wall bounding 
the fluid. However, when the system is a finite one, such a pressure would depend 
on the position and shape of this wall, and also, for a fixed curvature, the mere 
presence of a boundary would change the size of the system. In the case of a 
pseudosphere, a system with a boundary makes difficulties because its perimeter 
grows as fast as its area. The stress tensor approach, which has been used here, has 
the advantage of defining a bulk pressure which does not refer to any wall. But it 
raises a question: should the self part be included or not in the definition of the 
pressure? One might be tempted to follow the usual procedure of discarding self 
effects. However, as seen in Section 5, it is necessary to keep this self term for 
retrieving a finite pressure in the Coulomb case. Thus, the nonself and self parts 
of the pressure are somewhat entangled with each other. For a flat system in IR^ 
or 11?, the (properly regularized) self pressure vanishes, for a Yukawa system and 
thus in the Coulomb limit. However, the same pressure for a flat Coulomb system 
can be obtained by starting with a Yukawa system (plus background) in S^, S*^, or 
the pseudosphere, going to the Coulomb limit a ^ 0, and finally going to the limit 
i?— ^ooora^ooofa flat system; if this route is followed, it is mandatory to keep 
the self part of the pressure. 

The entanglement of the nonself and self parts of the pressure is especially appar- 
ent in the two-dimensional case. The result for a flat Coulomb system, Pcx = —nq'^/4 
can be obtained by starting with a Yukawa system on a sphere S"^, and taking the 
limits q; -^ and R ^ oo. These limits can be taken in one order or in the opposite 
one, or even for a fixed value of aR, giving the same final result for Pcx- However, 

12 



the separate contributions Pnonscif and Pscif do depend on the way the hmits are 
taken. 

In a curved space, the pressure from the stress tensor approach can also be re- 
trieved from the usual thermodynamic definition (1.2) (the negative of the derivative 
of the free energy with respect to the volume, or its two-dimensional analogs). If 
one deals from the start with a Coulomb fluid in S^ or S"^, there is some arbitrari- 
ness in the definition of the self-energy of a particle P], and only a "reasonable" 
choice allows to retrieve ||l|] the stress tensor pressure. In the case of a Yukawa fluid, 
the self-energy (3.9) or (4.10) is sufficiently well-defined for the calculation of the 
pressure. Starting with a Yukawa fluid and going afterwards to the Coulomb limit 
avoids the above mentioned arbitrariness. 

It should be noted that, in curved spaces, the stress tensor approach defines 
the pressure as the response to a change of the radius R or a; the thermodynamic 
approach does the same. However, as i? or a changes, the interaction potential 
changes, and therefore the derivative —dF/dV is not taken at constant interaction 
potential. This is a bit unusual! It has to be taken into account when using the 
virial expansion which seems to exist in the pseudosphere. 

The conclusion is that there is some arbitrariness in the definition of the pressure 
in a curved space. However, the stress tensor approach, in which there is some reason 
for including the self contribution, gives a pressure which seems acceptable. 
APPENDIX A. FORMULAS IN S3 

For deriving (3.3) from (3.2), one first note that, since the integrand in (3.2) 
depends only on the shape of the geodesic triangle formed by points (0,1,2), the 
integration can be performed on another pair of positions (0,2), rather than (1,2). 
By an integration by parts on 0, the term VoG'('?/^oi) " VoG'('?/^o2) can be replaced by 
—G{ipoi) Aq G{iIjo2)- Using the Helmholtz equation with a point source 

- Ao +«2] G(^o2) = -4:7rS^'\i^o2) (A.l) 

gives 

,2 



^2^2 r 

-* nonsclf 



247r 



47r / dV2G{iJi2)g{i^u) +2a' / dVodV2G{iJoi)Gii;o2)gi^i2] 
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(A2) 



Using the Dirac notation 



G{iJij) =< i\ 



An 



^J> 



A+a2' 
where i,j are positions, gives for the integral on in (A. 2) 

47r 



(A3) 



J dVoG{^,o)G{^o2) 



dVn < 11 



A+a2 



=< 11 



(47r) 



|2 >= -47r 



0x0 



dG{ij^2] 



A-n 



A+a2 



12 > 



(-A+a2)2i 5(^2) 

Thus, with simpler notations V and ip instead of V2 and ■?/'i2, (A. 2) becomes 

dG{ij) 



-Pnonself = ^ / C?V^ 



G(^)-a- 



9q; 



g{^) 



{A.A) 



(A5) 



Finally, since RG{iIj) as defined in (3.1) depends on a only through aR, a{dG/da) = 
d{RG)/dR and (A.5) gives (3.3). 

For deriving (3.6) from (3.4) and (3.5), first one performs an integration by parts 
which transforms (3.5) (where dV = An R^ siii^ ip dip) into 
.2 



Pi 



nq 

24^ 



l/)>l/)0 



-G(^) A G(V') + 3^2 [G{iP)Y 



dV 



dC 

G{ilj)4^AnRsm^ij 
dip 



ip=ip(, 
(A6) 

Since G obeys the Helmholz equation and ■?/' 7^ 0, in (A. 6) the Laplacian A can be 

replaced by a^. Furthermore, the integral remains convergent as V^o ^ and can be 

extended to the whole hypersphere. With G given by (3.1), this is an elementary 

integral on trigonometric or hyperbolic functions. Adding the savaXl-ipQ form of the 

last term of (A.6), one finds 



Pi 



6 



nq 



na^R 



cot TTo; 1 
+ 



+ 0{% 



for aR < 1 



(A.7) 



sin nuj Rjj RipQ 

Adding (3.4) and (A. 7) gives (3.6). A similar calculation can be done when aR > 1. 
For computing the self-energy (3.9) of a particle pictured as a spherical surfacic 
charge of radius tq = Rsinipo, in the small-V^o limit one can use the point-particle 
potential (3.1). As ipo ^^ 

uo cot nuo 



q^ q^ \ 1 

Cseif = yG(V^o) = y 



0{iPo 



for aR < 1 



(A^ 
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Up to the volume independent term g^/2ro (which becomes infinite in the ipo -^ 
hmit), (A. 8) does give (3.9). For aR > 1, uj must be replaced by iuj. 
APPENDIX B. FORMULAS IN Sa 

For finding the Green function (4.5), one notes that the Helmholtz equation in 5*2 
(similar to (A.l) with 27r5^'^^ instead oi AnS^^^) reduces to the Legendre equation [Q, 
with the variable cos^ and the parameters u = (l/2)[— 1 + (1— 40;^/?^)^/^], yU = O.The 
solution singular at 6* = and regular at 6* = vr is the Legendre function 

P.(-cosg)=F(-t.,z/+l;l; ^^^°' ) (5.1) 

Indeed, a.t 6 = tt the hypergeometric function F is regular, and as 6 ^ its behavior 
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p f n\ '^ sin uir 
RA— cos 0) = 



TT 



TT 

In sin — h 7 + ip{l + y) ^ — cot vtt 



0(1) (5.2) 



isQ such that (4.5) does behave like — In ^. 

For deriving the self-energy (4.10), one looks at the behaviour of (g^/2)G(6'o), 
considering tq = R9q as a small fixed particle radius. As 6*0 -^ 0, using (B.2) in 
(B.l) gives (4.10). 
APPENDIX C. GREEN FUNCTION IN A PSEUDOSPHERE 

The Helmholtz equation in a pseudosphere of "radius" a reduces to the Leg- 
endre equation p with the variable coshr and the parameters z/ = (1/2) [—1 + 
(1 + 4a^a^)^/^], /i = 0. The solution singular at r = and vanishing at infinity is 



Qu{cosh.T). For small r, it does behave [jTO[ like — Inr. 
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TThere is a minor misprint in ref. ||T0|. In eq.(15) p. 164 7 must be replaced by 27 for this 
equation to agree with eq.(12) p. 110. 
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